// distance.sci
  function dist = distance(x,y,z)
    dist = sqrt(x*x + y*y + z*z)
  endfunction